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Special case calculations are presented, which can be used to put restrictions on the general form 
of heat kernel coefficients for transmittal boundary conditions and for generalized bag boundary 
conditions. 
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The heat equation asymptotics plays a prominent role in mathematics p7[ | and theoretical physics 
^ IiUjIiJiQ- For manifolds without boundaries efficient schemes have been developed and the calculation 
of the volume part is nowadays nearly automatic |26 |,^^- When the manifold has a boundary, suitable boundary 
conditions have to be imposed and the heat equation asymptotics receives additional contributions depending on 
• the boundary conditions considered . In recent years a conglomerate of methods has been proven to be very 

effective in the determination of this asymptotics for Laplace type operators. The conglomerate of methods consists 
^ of functorial techniques, index theorems and special case calculations. All ingredients are very agil in that they can be 
adopted to the specific boundary condition considered. A characteristic feature is that on its own none of the methods 
I is able to determine the coefficients fully, but instead each determines only part of the complete answer. As a rule, 
the information obtained by each ingredient overlaps with the information obtained from the others and this provides 
T— I . crucial cross checks for the calculation. Roughly one could say, that special case calculations are "responsible" for 
^ ' the group of extrinsic curvature terms, the index theorem allows to extract information on the gauge connection 
\l ' terms and functorial techniques, most prominently conformal transformation techniques [0, determine whatever is 
left over. Employing this procedure, the heat equation asymptotics for the classic as well as for "exotic" boundary 
conditions has been determined. A summary is provided in the conference contribution As further approaches 
^ . let us mention (3^||,|^,|l§. 
^— I ■ The special case calculations that have been employed so far, consist of the following examples. (1) Dirichlet 
and Robin boundary conditions: scalar fields on the generalized cone / x M, I — [0, 1], A/" a Riemannian manifold 
|7|,p| |3^j2^ . (2) Mixed boundary conditions: spinor fields and p- forms on the generalized cone 

(3) Oblique 

> boundary conditions: again the generalized cone, in addition x T^^^ with the two-ball B and the {D — 2)- 
dimensional torus T^^^ |^^. (4) Spectral boundary conditions: generalized cone plus B^ x J\f |2l| , [l9| , p9| . 
I In the present contribution we will present two further special case calculations, one on a different manifold and one 
i on the ball with different boundary conditions. First we provide a calculation on a manifold obtained when glueing 
^ ' together a hemisphere and a ball. This example determines crucial information for transmittal boundary conditions 
• i-H pH . Second we consider the Dirac operator on the ball with boundary conditions as they occur in gauge theories 
in Euclidean bags. These boundary conditions involve an angle 9, which is a substitute for introducing small quark 
H , masses to drive the breaking of chiral symmetry [ ^3p^ . An analysis of the associated heat equation asymptotics has 
- - . barely started ||,|8|. 



II. HEMISPHERE BALL GLUED TOGETHER 



The first example we present is of relevance in the context of transmittal boundary conditions. To define transmittal 
boundary conditions, see e.g. consider the d-dimensional manifold M = M.^ Us , which is the union of 
two compact manifolds M.^ along their common boundary S. Let be Laplace type operators on M.^ written 
invariantly as 
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The operator D = {D~^,D ) acts on a pair = {(j>'^,4> ). The transmittal boundary condition is defined by the 
operator 

Bu^ = - rls} © {(V+ + 0+)|s + (V„_r )|s + C/0+Ie}, (1) 

where is the respective exterior normal derivative on to the boundary S. 

Let us consider a specific example of the above situation and we start with the two-dimensional case. 

Let — be the unit hemisphere and A^^ = the two-ball. On the hemisphere we consider the "massive" 
Laplacian 

D+ = -Ah2 + J. 

Adding the curvature term (l/8)i? = (1/4) has the advantage that the eigenvalues of Z?+ are complete squares. In 
detail, the eigenfunctions are 

<l^herms{0, if) = iV/.e^'^F^^i^/l^ (cos 0), m € 7L, 

with the associated Legendre functions and Nh a normalisation constant. The values of A have to be fixed by 

boundary conditions. 

On the ball we simply take 

D- = -Ab2, 

with the well known eigenfunctions 

<pd^sc{r, if) = Nde"'"^Jim\ (Ar), m e TL. 

Here, Jv{x) are the Bessel functions and Nd is a normalisation constant. 

Next we glue the hemisphere and the disc together along their boundary, which is a circle. As a specific example 
of the transmittal boundary conditions (|l|) we impose that the eigenfunctions of and D~ as well as their normal 
derivatives agree along the circle, this is we choose [/ = 0. 

Matching the eigenfunctions eliminates one of the normalization constants, 

m| 



so 



Nd = A^Pj7l2(0), Nh = iVJ|,„|(A). 

Next we match the normal derivatives. The natural normal derivatives in the example are the exterior normal 
derivatives [d/dr) and {d/d9). In these the condition reads 

oO or 



and this gives the implicit eigenvalue equation 

= \Jl^^iX)p-%{0) - J|„|(A)^P,-i™l,(x) U=o. (2) 
In the following we will simplify the notation to 

Once this implicit eigenvalue equation (|^) is known, the procedure developed in ^ can be used in order to determine 
the heat equation asymptotics for the specific problem at hand. Starting point of the procedure is the contour integral 
representation 



(3) 
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where the contour 7 encloses counterclockwise all real positive zeroes of equation (||) and where a small mass (as 
a infrared regulator put to zero later) has been introduced. 

The next step is to deform the contour to the imaginary axis, obtaining 



00 ^ 



(4) 



where P^{x) = P^^^_-^{x), this implies Pii}^i/2i^) — been used ]3^ . 

As explained in great detail in [Q, the heat equation asymptotics is encoded in the asymptotic behavior of the 
implicit eigenvalue equation (^) . For the Bessel functions everything needed is provided by Olvers uniform asymptotic 
expansions. For m — > 00 and z — k/m fixed, they read 



Imimz) 



V2^(l + z2)i/4 



1 + E 



Uk{t) 



fc=l 



1 £'""(1 + z2)l/4 



1+E 



Vk{t) 



(5) 
(6) 



with t = l/vT+z^ and 77 — vl + z^ + ln[z/(l + \/l + z^)]. The polynomials are defined recursively through 

t 
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Uk+iit) = -t'il - t'Kit) + - / dr (1 - 5T')ukiT), 



Vk{t) = Uk{t)+t{t^ - 1) 



-Uk-l{t) + tu'f,_^(t) 



starting with Uf^{t) = 1. The first few coefficients are listed in Q, higher coefficients are immediate to obtain by the 
above recursions using a simple computer program. 

The relevant information for the Legendre functions is p3] 



r(3/4 + (m + A)/2)r(3/4 +{m- A)/2) ' 



r(l/4 + (m + A)/2)r(l/4 + (m - A)/2) ' 
Further expansion of the F-functions is achieved using [3^ ] 

B2k 

2fc(2fc- l)z2fe-i' 



(7) 
(8) 
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with the Bernoulli numbers Bik- We use this in the slightly different form 



r(u + a) ^ exp [ ulnu — u + (a — 1/2) Inw + ln\/27r + 



7/ (a) 



i=\ 



respectively 



1 / 00 

— ~ exp I — wluM + u — (a — 1/2) In u — ln\/27r ) 1 + 7 

F(u + a) V ^ \ ^ 



7i(a) 



with the 7; (a) easily determined by an algebraic computer program. 

These expansions are used for to 7^ after substituting k —> km in the integral (Q) ; a suitable choice turns out to be 
u = to(1 + ik)/2. As known from the ball, to = needs special treatment and a simple large argument consideration 
is sufficient. 
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Without going into the details of this calculation, as to be expected the sum over m produces Riemann zeta functions 
and a typical expression for the asymptotic contribution is 

oo 

TT J ok 



where, e.g.. 



24(l + fc2)3/2 ' ' 32(l + /c2)3- 

The fc-integrals are standard and residues and values of C(s) are easily determined. As a pleasant feature let us 
mention that the whole calculation can be and has been fully automised with a simple algebraic computer program. 

The transition from two dimensions to arbitrary dimension d is surprisingly easy once the suitable organising 
functions are known. On ball and sphere settings the best organisation is obtained via the use of the Barnes zeta 
function 

CMS, b) = f: 77 r = flCY'.^)ii + br'- 

t,i^ + ^^i + - + m,Y d-1 J' 

In the present context this is seen as follows. On the hemisphere, eigenfunctions can be represented as Gegenbauer 
polynomials ||3^ . The relevant dependence on the normal coordinate 9 is 

where again a coupling 

has been included. The radial dependence on the ball is 

■/i+(d-2)/2(Ar) 
r(rf-2)/2 

such that 

d Ji+(d-2)/2(Ar) d~2 

^(rf-2)/2 - AJ,+(d_2)/2(Ar) —Ji+{d-2)/2[>^r). 

The angular tangential dependences of both solutions agree. Eq. (||) suggests, that the most convenient choice to 
proceed is the implicit eigenvalue equation 

- AJ,V(,_„/2(A)Pr'rr'^/^0) - Jz+(.-2)/2(A)^P,-i7;^^)/^(0). (10) 
The index m is replaced by ^ + (d — 2)/2 and the degeneracy 



d{l) = (2/ + d-2)- 



ll{d-2)l 



of the spherical harmonics has to be taken into account. Thus final answers can be found from the two dimensional 
result once the Riemann zeta function is replaced by a sum of the Barnes zeta function , 



1=0 



E <0 (l + ^j = Cb(2s, {d - 2)/2) + C6(2.s, d/2). 



Note, however, that eq. (10) does not correspond to matching the normal derivatives at the boundary, but instead 
they are assumed to have a jump described by[/=(rf— 2)/2, 
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dr 



<pdi. 



r=l + 



89 



(f>he 



d-2 



e=Tr/2 



-(f>d'is 



In summary, the meromorphic structure of the zeta function for the d-dimensional problem is completely clear and 
can be used to calculate the leading heat kernel coefficients. Subtracting the volume contributions of and , 
the "boundary part" of some leading coefficients is 



ai = 0, 



02 



as 



04 



(4 - d)2-'' 
3r(d/2) ' 
(d-5)(d-3)0F2-'' 



64r(d/2) 



180r(d/2) 



250 2839 , 191 61 

\ d d^ H d-^ 

7 42 7 21 



These results can and have been used to put restrictions on the general form of the coefficients for transmittal boundary 
conditions pl[ |. To exemplify the procedure consider the leading coefficients. Let be the extrinsic curvature on E 
as induced from Ai^. Invariance theory shows that for / a localizing test function 

= Jdy Ci/Tr (1), 
s 

= (4^)-'*/2i J dy Tr (ei/(if+ + R-)! + e2(/+ + + + e^fu) . 

s 

For the hemisphere ball example, we have K'^ = 0, = d — 1, U — {d — 2)/2, and derive the equations 

ci -0, 

(d-l)ei + — — 63 = A-d, 



and so ei = 2, £3 = —6. This information, as well as the remaining one £2 = 0, is easily obtained by different means, 
but for the higher coefficients the special case input is crucial. 



III. GENERALIZED EUCLIDEAN BAG BOUNDARY CONDITIONS 



In theories of Euclidean bags, chiral symmetry breaking is triggered by imposing the boundary condition 

= n_v^ \dM 

■.= l[l + ^e<>^''T'T"^)\eM (11) 

on the spinor field -ijj. For 9 — this is a well studied problem which can be understood as a mixed boundary problem 
of an associated Laplace type operator D. Let ej be a d-bein system, Tjki the Christoffel symbols relative to the 
orthonormal frame and F-' the F- matrices projected along this frame. Let P be the Dirac operator on M, 

P = -iV^Vj, 

where Vj — Cj + tOj is the covariant derivative with the spin connection 

ujj = --Fj-feiF'^F'. 

Consider the associated second order problem for D = P^ with domain 

domain^) = {V G C^{V) : Il^iP \aM ® n_ {-ij'Vj) V-jg^ = 0}. 
For = 0, this leads to the mixed problem 



5 



(12) 



where 



and 



n+ = -(i- zr^r" 

+ 2 ^ 



In order to derive eq. (p^), it is crucial that [n-jF"] = 0, see a relation that does not hold for 0^0. For the 
case = the special case calculation on the ball has been performed in jl^,^ and subsequently used to restrict 
the general form of the coefficients for mixed boundary conditions pl| . It is the aim of this section, to generalize the 
special case calculation to ^ and to use the results in order to determine the oi-coefficicnt for a general manifold. 
The structure of the higher coefficients is considerably more difficult due to the non-commuting F-matrices involved 
in the invariance theory. However, work along the lines shown here is in progress ]2^ . 

We start the analysis by solving the eigenvalue problem on the d-dimensional ball. The eigenspinors of the Dirac 
operator on the ball have the form |jl^ 

r('i-2)/2Uj„+(,_2)/2(A:r)4")(f]) 



with n e INo. Herc^^z£"^ (f7) are the spinor modes on the sphere, see [|l^, and C is a normalisation constant. 
The condition (y), 

is easily applied to (|l^) and (0), and the implicit eigenvalue equations read 

Jn+d/2{k) T e'^ Jn+d/2-l{k) = 



for V'i^^ , and 



Jn+d/2{k) ± e ^ Jn+d/2-l{k) = 



for Suitably combined, we write these as 

Jl+d/2-i{k)-e-''Jl+„^{k)=0, (15) 
Jl+d/2-i{k)-e''jl^„^{k)=0. (16) 

Clearly it is sufficient to deal with cq. (|l5|), the contributions from (|l^) follow by replacing —0. 

For convenience we introduce v = n -\- [d — 2)/2. To analyse the zeta function we start again with the contour 
integral representation 

C(.)=^d(.) j ^fc-^^Ain(j2(fc)_e-2«j2^^(fc)), 

7 

where d{i') is the degeneracy associated with the implicit eigenvalue equation ( [l5| ) and 7 again encloses counterclock- 
wise the positive real zeroes of (p^). (For convenience, we omit writing the infrared cutoff Af^. The procedure in 
principle is as in the previous section.) Denoting by ds the dimension of spinor space, dg = 2^*/^, and taking into 
account the sphere eigenspinor degeneracies, one finds 
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Proceeding in the manner described, we shift the contour to the imaginary axis to find 

oo 

C(^) = ^ E '^(-) / fc-^^^ In im + e-'Oll^,{k)) . (17) 



To simplify the analysis of the asymptotic behavior of the integrand it is convenient to rewrite ( p^ in terms of Bessel 
functions involving only one index. For this purpose we use p2| 



which allows to write 



oo 



^ In [e-'ll\zv) + ll{zv) - h^' Uzv)ll{zv)^ , 

when irrelevant factors in the logarithm are neglected. 

The uniform asymptotics is completely determined by eqs. and (^. Here we will concentrate only on the terms 
contributing to the heat equation coefficients oq and ai. It is easy to see, that the relevant pieces from the argument 
of the logarithm are 



(...) . ^ - 2(1 - 1) |i + ^ r - ij I + i-. (18) 

The contribution to C(s) resulting from the first line has already been dealt with in the calculation for the case 9 = 
and it reads ||2^ 

4 ris-i/2) ^ ds r(s + i/2) 

4V7r i(s + l) 4:^/^T l(s + l) 

The same contribution comes from (|T6|). The second line in ( |l8| ) contributes the ^-dependent piece 

C,is) = ^ E di.)u-'^ J dz z--^ In (l + [e^^ - l]) . (19) 



Using standard Taylor series expansion, one first finds 

|l„ (l + if r - 1]) -g-'ft"-"^)"' ,!,, ^ 

This allows the z-integrals to be done and an intermediate answer is 

^g(,s) = ^Cg(2-s, d/2 - 1) y (-l)'(tanhg)'+i + ^,^{^^ . 

4r(s)^^^ ' ' ^ ' r((/ + 3)/2) 



Simplifications occur if C,_e from (|T^ is added. 



[ 

^C/3(2s,d/2- l)2i^i(l,s + l,2,tanh 
^C6(2s,d/2-l)(cosh2^0-l). 



Ceis) + C-o{s) = ^C6(2s, d/2 - 1) ^(tanhfl) ^ 

2 /l\ 
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As expected, ior 6 = the new term Ce{s) + C-e{s) vanishes. From here, the residues Res C{d/2) and Res C{{d— l)/2) 
are easily determined. The normalization coefficient uq is confirmed, for ai we find 

ai = (47r)-('^-i)/2i4|^'-'| (cosh'^-i0- 1 
with 15''^^ I the volume of the {d — l)-sphere. On an arbitrary manifold M. the coefficient is written in the form 



dM 

and we read off 



s = -( cosh'*"^ e - 1 

4 V 

As expected, the coefficient S depends on 9 as well as on the dimension d. The dimension dependence is a result of 
the fact that the number of F-matrices depends on the dimension d. 



IV. CONCLUSIONS 



This contribution provides further evidence that special case calculations play an important role in the determination 
of the boundary contributions to the heat equation asymptotics. 

Whereas previous applications involved only Bessel functions as special case solutions, in our first calculation we 
have provided here an example where the associated Legendre functions are of relevance. The example involved was 
the union of a hemisphere and a ball of radius one, as it turned out to be useful in the analysis of transmittal boundary 
conditions |3^] . The choice of a hemisphere simplifies the calculations considerably, because the asymptotic behavior 
of the Legendre functions needed is relatively simple, see eqs. and (^. If instead of the hemisphere a spherical 
cap is considered, further contributions will arise because the boundary is not geodesically complete any more and so 
Kab ^ 0. Presumably, also for this more general situation an analysis should be possible, for example along the lines 
of 1^, where the needed asymptotic behavior of the Legendre functions has already been determined. 

The nice feature of using the zeta function for the analysis of the heat kernel asymptotics is that it also allows for 
the consideration of other spectral functions like the determinant and the Casimir energy. So as a side effect of the 
present calculation an investigation of the influence of edges on these quantities can be envisaged. 

In the second application presented, we have provided the first steps into the thorough analysis of the boundary 
condition (pi]). The explicit determination of the heat kernel coefficients is in its infancy. A special case calculation 
is sufficient to obtain the full ai-coefficient, a result which is already new. For the higher coefficients, techniques 
particularly adapted to the Dirac operators are needed ||2^ and further results will be presented in p^ . 
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